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Abstract 

We propose the minimal, lepton-number conserving, SU(3) c xSU(2)ixU(l)y gauge- 
singlet, or phantom, extension of the Standard Model. The extension is natural in the 
sense that all couplings are of 0(1) or forbidden due to a phantom sector global U{\)d 
symmetry, and basically imitates the standard Majorana see-saw mechanism. Sponta- 
neous breaking of the U(l)n symmetry triggers consistent electroweak gauge symmetry 
breaking only if it occurs at a scale compatible with small Dirac neutrino masses and 
baryogenesis through Dirac leptogenesis. Dirac leptogenesis proceeds through the usual 
out-of-equilibrium decay scenario, leading to left and right-handed neutrino asymme- 
tries that do not fully equilibrate after they are produced. The model contains two 
physical Higgs bosons and a massless Goldstone boson. The existence of the Goldstone 
boson suppresses the Higgs to bb branching ratio and instead the Higgs bosons will 
mainly decay to invisible Goldstone and/or to visible vector boson pairs. In a repre- 
sentative scenario, we estimate that with 30 fb _1 integrated luminosity, the LHC could 
discover this invisibly decaying Higgs, with mass ~120 GeV. At the same time a signif- 
icantly heavier, partner Higgs boson with mass ~210 GeV could be found through its 
vector boson decays. Electroweak constraints as well as astrophysical and cosmological 
implications are analysed and discussed. 



1 Introduction 

The Standard Model (SM) has just two openings where renormalisable operators can be 
added which couple SU(3) c xSU(2) L xU(l)y gauge singlet fields to SM fields. One place is 
the super-renormalisable Higgs mass term £Q, the other place is the lepton-Higgs Yukawa 
interaction LH. What would happen if we filled in these gaps? 

There is no physical evidence, as yet, to suggest that B — L (baryon - lepton number) is 
not a good symmetry of nature. The SM preserves B — L, so we will choose to extend the 
SM in a B — L preserving way. However, overwhelming evidence supporting small, non-zero 
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neutrino masses does exist |2j. We are therefore led to build a model with Dirac masses 
for the neutrinos and see if it is possible to create the observed baryon asymmetry of the 
Universe within this set up. Ideally, we should also strive to build a natural model, both in 
the 't Hooft sense and the aesthetic sense. In particular, Yukawa couplings should be either 
(9(1) or strictly forbidden. 

Following this approach, we augment the SM with two SU(3) c xSU(2) L xU(l)y gauge 
singlet fields, a complex scalar $ and a Weyl fermion sr. These fields will provide the link 
between the SM and a phantom gauge singlet sector. 

-C hnk = (k^-Hsr + H.c.) - rj H^H $ , (1) 

where h v and r\ are dimensionless couplings of 0(1) in line with our naturalness criterion. The 
field H (or H — ia 2 H*) is the standard SU{2) L Higgs doublet 1 responsible for spontaneous 
electroweak, SU(2)l x U(1)y, gauge symmetry breaking. Note that sr must carry lepton 
number L = 1. 

In this form, the model is incomplete because neutrinos would have large, electroweak- 
scale masses. However, we have so far ignored the purely phantom, gauge singlet sector of 
the model. Here we add a phantom right-handed neutrino ur, and sl the partner of the 
Weyl field sr. These fields will also carry lepton number. The fermionic part of the phantom 
sector therefore contains 

- £ v = h v $sZv R + MsZsr + H.c. (2) 

where h p is a general complex Yukawa coupling of 0(1). Other possible lepton number 
conserving terms, 

T L Ev R + M'sZur + $sZsr + H.c. + ... , (3) 

are forbidden when imposing a global U(1)b symmetry, under which only the fields, 

VR-^e ta u R , $^ e -' ia $, (4) 

transform non-trivially. This choice for the phantom sector is purely motived by the need for 
the simplest model leading to small, Dirac neutrino masses. A crucial point to notice here is 
that the spontaneous symmetry breaking of the global U(1)-d will trigger consistent electroweak 
gauge symmetry breaking through the last term in Aink provided that ($) = o ~ v , with v 
being the vacuum expectation value (vev) of H. 

The model 

£■ = £sm + ^ link ' (5) 

can be trivially embedded into Grand Unified Theories (GUTs) and can be supersym- 
metrized. Part of the model was first presented in the literature by Roncadelli and Wyler 31 . 
who were motivated by the need for a model with naturally small Dirac neutrino masses jlj. 

l ln our notation is H T = (H+,H°) and H T = (H°* , -H~). 
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Figure 1: The diagram responsible for light Dirac neutrino masses in the model 

Notice that its structure is different (and much simpler) than the ones exploited recently by 
EH E] - though the latter are supersymmetric. 

We will confine our discussion to a three generation neutrino model. We will therefore 
add three generations of the sr, sl and vr. For simplicity we will consider just one copy 
of the complex scalar $. This is the Minimal Phantom Sector of the SM consistent with 
naturally small Dirac neutrino masses, and as we shall see shortly, provides an explanation 
for the baryon asymmetry of the Universe and has Higgs phenomenology strikingly different 
to that of the SM. 

2 Neutrino Masses 

In this section we briefly repeat the main points of [B]. Notice that we are free to work 
in the basis where the Dirac mass matrix M in (J2J) is diagonal. After spontaneous £/(1)d 
symmetry breaking, the Lagrangian (jSJ) results in Dirac-neutrino effective mass terms of the 
form, m y z/ R + s' L niN s' R , where up to terms (9(M~ 2 ), we obtain 

m„ = — m M -1 m p , rriN = M , (6) 
with m = h u v and m p = h p a being 3x3 matrices and with neutrino mass eigenstates 
v' L = z/ L - mM x s L , s' L = M 1 v\, + s L 
*4 = - mpM^SR , s R = M^mpZ/R + s R . (7) 

Bold face letters denote 3x3 matrices and column vectors. From © we obtain a typical 
seesaw spectrum with light and heavy Dirac neutrinos with masses and M, respectively. 

The physical neutrino masses will then result from the final rotation = A'l'm^B where 
A, B are unitary matrices. In a basis where the charged lepton Yukawa couplings, h e are 
diagonal, the matrix A will just be the usual PMNS matrix jS], measured by neutrino flavour 
oscillation experiments. Minimal mass matrices for Dirac neutrinos have been recently 
classified in jHj and can be exploited to shed light onto the connection between CP-violation 
and Dirac leptogenesis that follows in the next section. The reader has possibly already 
realized that the model in © is a simple realization of the Froggatt-Nielsen ^Oj mechanism, 
usually invoked to generate large hierarchies in quark masses. In this model neutrino masses 
are generated by the Feynman diagram in Fig. ^ and are given by 

il% = 9ij e a ^ , (8) 
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with e = cr/M and Oj + bj = 1, corresponding to the difference between the [/(l)D-charges 
of the left and the right handed neutrino. The matrix is a general matrix, in our case 
a product of the matrices h u v and h p . In contrast with the quark case where e naturally 
explains the large hierarchy in quark masses, here the smallness of e explains the relative 
smallness of neutrino masses. Assuming that the Yukawa couplings are perturbative and of 
order one, with v = 175 GeV and m v = 0.05 eV we find the ratio e = cr/M, needs to be 

e ~ 3 x 10" 13 . (9) 

In the next section we will see that it is possible to achieve this in a model consistent with 
our naturalness criterion and successful baryogenesis jllj . Significantly, in this model there 
is no neutrinoless double beta decay, however, we will later examine the consequences of the 
phantom sector for Higgs searches at the LHC. 

3 Baryogenesis 
3.1 Dirac Leptogenesis 

Although B — L is preserved exactly in this model, we will see that baryogenesis through 
(Dirac) leptogenesis is still possible [TT]. Just as in the SM jHj, in this model the combina- 
tion B + L is anomalous and at low temperatures B+L violation proceeds through tunnelling 
and is un-observably small. At higher temperatures, close to and above the critical temper- 
ature for electroweak symmetry breaking, T c > 150 GeV, thermal fluctuations allow field 
configurations to pass over the 'sphaleron barrier', leading to rapid B + L violation [To] . 

It is important to note that the rapid B + L violating processes do not directly affect 
right-handed gauge singlet particles. Large Yukawa couplings between the SM quarks and 
charged leptons will however tend to equilibrate asymmetries in the left and right sectors 
of the model, depleting an overall 'right-handed' B + L as an overall 'left-handed' B + L is 
depleted via 'sphaleron effects'. 

The crucial idea behind Dirac leptogenesis (or Dirac neutrinogenesis) is that the small 
effective Yukawa couplings between the SM Higgs and the left and right handed neutrinos 
could prevent asymmetries in the neutrino sector of the model from equilibrating. Therefore, 
even in a model conserving total lepton number, a left-handed B — L asymmetry could be 
produced at the same time as an opposite, right-handed B — L VR asymmetry. The left-handed 
B — L asymmetry would then lead to an overall baryon asymmetry just as in Majorana 
leptogenesis. Clearly, for this to work the effective neutrino Yukawa couplings must be small 
enough to keep the left and right lepton asymmetries from equilibrating until (at least) after 
the electroweak phase transition, when the sphaleron processes linking the baryon and lepton 
asymmetries would have dropped out of thermal equilibrium. This mechanism even works 
when the initial, overall B = L = 0. It is especially interesting to note that this mechanism 
links the baryon asymmetry directly to the smallness of the Dirac neutrino masses. 

At temperatures above T c , when sphaleron and other SM processes can maintain most 
SM species in thermal equilibrium it is possible to derive relations amongst the chemical 
potentials of the various particle species [T^]. Since we demand the right handed neutrinos 
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Figure 2: Feynman diagrams for the decay of the heavy gauge singlet Si into a v Rk and a 
$. The CP-asymmetry in this decay is due to the interference between the tree-level diagram 
(a) and the 1-loop self-energy diagram (b) . 



be out of thermal equilibrium we can ignore their contribution for the moment, leading to 
the usual relation between baryon and lepton number used in Majorana leptogenesis 

Y B = ^(Y B -Y LsM ), (10) 

where Yx refers to the 'asymmetry in X' to entropy ratio and LgM refers to the lepton number 
held in SM particles (not including the right handed neutrinos). If we now consider the case 
where, initially, the total B — L was zero, we have the relation 

Y B - Y LsM - Y Lvr = 0, 

which in conjunction with equation fTT)|) leads to 

28 

Y B = -t^Y l „ r , 

showing that just over a quarter of the right handed neutrino asymmetry is 
a baryon asymmetry. 



(11) 
(12) 

converted into 



3.2 CP-violation 

As is well known, in order to generate a particle-antiparticle asymmetry in the early universe, 
the three Sakharov criteria must be fulfilled ^Tj\. Particularly relevant to this model are the 
requirements for a departure from thermal equilibrium and CP-violation. 

In analogy with conventional leptogenesis, in this model, CP-violation and a departure 
from thermal equilibrium could arise during the decays of the heavy Dirac S = sl + sr 
particles. In particular, CP-violation would originate through the interference between a 
tree level decay and a 1-loop self-energy diagram as shown in Fig. |21 If we define a 'right' 
CP-asymmetry as 

Ek (r(5i - i/fl fc $) - r(5, - 9 Rk ®*)) 

5 Ri = ^ jT (s i ^u Rj *) + '£ i nS i -+L l H) ' (13) 
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where T(Si — > VRk$) is the rate of the decay process Si — > VRk& etc. L% and H repre- 
sent lepton and SM Higgs SU(2)i doublets respectively. In addition, unitarity and CPT 
conservation provide the following useful relation 

Y ?{Si - v Rj $) + Yl ros -^l 1 h) = j2 r os - vr? $*) + r os Li ' Ft ) ' ( 14 ) 

j i r v 

which leads to the relation 5 Li = —Sri, where 5n is defined in analogy with equation (JT3j) . 
We also define right and left branching ratios as 

5/ " = E.r^-^^ + ^r^-L^) ' (15) 

and B Li = 2- B Ri . 

In the limit that the Mj are hierarchical, the lepton asymmetry is generated via the decay 
of the lightest S and S. In this case, the CP-asymmetry (|T3j) is given by 

1 Mi Im (hpht)y (ht h u ) n 
^^8^1. W 3 (h p ht) 11 + (hth,) n ■ (16) 

It should be noted that the structure of the amplitudes leading to CP-violation in this sce- 
nario are similar to the self-energy contribution to the CP- asymmetry in Majorana leptoge- 
nesis. Therefore, in the limit of quasi-degenerate masses for the Si, a resonant enhancement 
of 5r should be possible [TBI CHI- We will however leave this case to be considered elsewhere. 
Using equation (jHJ), (fTH|) can be rewritten 



1 Mi Im (K m t h v)n 

5™- ~ 



"Ti v a (hphj)n + (hth,) n 



(17) 



Following the approach of j2Q], the most general h„ and h p can be parameterised in the 
following way 

^ = ^ AD ^ WE Vici> ( 18 ) 
h P = ^ D v^ XtD ^ Bt ' ^ 

where W and X are general 3x3 matrices satisfying the condition WX^ = 1 and = 
Z for the diagonal matrix Z. This allows 5r\ to be written as 



1 Mi Im 

Ori ~ 



(20) 



8tt v a (Xt m„ X) n + ( Wt m u W) n ' 

It is now straightforward to show, in analogy with [21], that \5m\ is bounded from above by 

1 Mi 
16n v a 



\$Ri\ ~ TcZ ~7Z ( m ^3 — rn v 1 ) ■ (21) 
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Again, we should stress that this bound could be violated grossly when the Mi are nearly 
degenerate, and that this caveat should be noted when this bound is used later. As we seek 
the most minimal model, without additional flavour symmetries in the phantom sector, we 
will not consider this resonant scenario here. 



3. 3 Out of equilibrium decays 

In the early Universe, if the expansion rate is faster than the interaction rate of a particle 
species then this species can become decoupled from the thermal bath. This statement can 
be quantified by considering the ratio of the particle interaction rate T(T), to the Hubble 
expansion rate H{T). If this ratio is less than 1, then the species evolves out of thermal 
equilibrium [2~2"j . 

If we require successful Dirac leptogenesis, the distributions of left and right handed 
neutrinos should be prevented from coming into thermal equilibrium from the moment a 
B — L UR asymmetry is created until after the electroweak phase transition, when the rate for 
B + L violating processes will be much smaller than the expansion rate of the universe. 

Processes leading to left-right equilibration include LH <-> $ vr mediated by the s- 
channel exchange of an Si, L vp <-» H $ and L$ <-» vp H mediated by the t-channel exchange 
of an S^ Approximately, at high temperatures these processes have a rate 

r^(T)~^^T 5 , (22) 

which should be compared to the Hubble parameter in the relevant radiation dominated era 



H{T) = J 8 ^ I- , (23) 
v ; V 90 M P ' K J 

where g* ~ 114 is the effective number of relativistic degrees of freedom in the SM plus 3 Vr 
and 1 complex $, and Mp = 1.2 x 10 19 GeV is the Planck mass. The strongest constraint 
will come from the highest temperatures when T ~ Mi, i.e. those at which the asymmetry 
is generated; 

h u \ 2 \h p \ 2 < J_,/8^7 (24) 



Mi ~ M P V 90 

To more accurately consider this constraint we need to solve the appropriate Boltzmann 
equations (discussed next). The dominant contribution to left-right equilibration will come 
from the inverse decay and subsequent decay of a real Si or Si. 

In this model, a L VR asymmetry can be generated via the standard out-of-equilibrium 
decay scenario, in analogy with various GUT baryogenesis scenarios and Majorana leptoge- 
nesis [T2"l IT3*] . In the following, we will consider the asymmetry to be generated solely during 
the decays of the lightest S. Pre-existing asymmetries from, for example, the decays of the 
heavier Ss will be treated as initial conditions. 

We will assume at T > Mi the abundance of Si and Si is thermal (we will relax this 
assumption later), and the number density of Si, ns 1 — n§ 1 ~ n 7 , where n 7 is the number 
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density of photons. As the Universe expands and cools the number density of the Si (Si) 
must rapidly decrease if they are to remain in thermal equilibrium below T ~ M\. If the 
interactions allowing this (primarily S\ (Si) decays) are slow compared to the expansion rate 
of the Universe then the Si and Si abundances will depart from their thermal equilibrium 
values. When the 5*1 (Si) eventually decay, the rates of back- reactions such as inverse 
decays will be suppressed by the relatively low temperature T <C Mi and the resulting L UR 
asymmetry will be [22] 



Y, 



We can define a parameter K such that 



K 



T(Si^v R §)+T(Si^LH) 
H(T = Mi) 



Sr n Sl 

Qjf Thry 



:h P ht)n 



6r 
9* 



(h+h 



(25) 



i/Jii 



Mr 



16nMi V 8tt 3 ^ 



90 



(26) 



where K <C 1 signifies that 5*1 is completely out of thermal equilibrium at T = M\\ the 'drift 
and decay' limit. We can re-cast this constraint in terms of 'effective neutrino masses' [2B] 
to make the connection with light neutrino data more transparent. Defining the effective 
neutrino mass as 



in 



> p ht)ii + (h£h„)u 



v a 
Wi 



K v a 



Mr 



IQtx 8n 3 g* 



90 



(27) 



we see that K < 1 is satisfied for m < m* where 



16vr /8 



v a 



Mr 



90 



(2i 



The connection between m and the physical light neutrino masses is clearly model dependent, 
and most applicable when (h p hj,)n ~ (h^h v )n. 

Finally, we can also introduce an efficiency parameter k, such that the L VR yield is given 

by 

(29) 



Yr 



9* 



For an initially thermal population of Si and when K <^ 1 the efficiency k ~ 1. 

If K > 1, Dirac leptogenesis can still be successful if the CP-asymmetry 8r is large 
enough. However, the simple estimate of the lepton asymmetry, equation (|23|) will no longer 
be valid, and the Boltzmann equations (BEs) should be solved. There are 4 coupled Boltz- 
mann equations relevant to this scenario, two for the Si total abundance and asymmetry, one 
for the asymmetry in and one for the asymmetry in L UR . Lepton number conservation 
means that only three of the asymmetry BEs are independent, we choose to eliminate the 
equation for L UR . The derivation of the BEs has been extensively covered in the literature, 
see for example [2312111221, so we will just write down the set of simplified equations for this 
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scenario. The BEs read 
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(30) 



where rjxs = (n s + n§)/n^, rj AS = (n s - n§)/n^, z = M 1 /T and 



7 " 



(2tt) 3 £ 



Si 



r(Si -> + 1(5! -> LH) 



-KAz) 



(31) 



where = 2 is the number of internal degrees of freedom of Si, Es 1 = \/p 2 + Mf and 
Ki(z) is a 1st order modified Bessel function. Decays and inverse decays of Si and Si are 
included through terms proportional to T D1 . Notice that these terms also include the most 
important CP-violating 2^2 scattering contribution coming from the subtraction of real 
intermediate states from the process LH <-> i/r<3> |24j . This subtraction is necessary to ensure 
unitarity and CPT are respected by avoiding double counting of processes in the BEs. 

T w parameterises the 'wash-out' due to processes of higher order in the Yukawa couplings 
which will tend to equilibrate the L and vr asymmetries, after the above subtraction of 
possible real intermediate states has been carried out. These process will be predominantly 
mediated by the off-shell exchange of Si^.3 and will therefore be highly model dependent 
[221 122] • For small Yukawa couplings where \h p \ 2 ~ \h^\ < 1, the processes LVr <-> H $ and 
L$ «-+ vr H will be negligible compared to decays and inverse decays. The contribution from 
the off-shell process LH «-> z/r$ is bounded from above by T D1 for the region around z ~ 1, 
if the Si have a reasonably large hierarchy in mass. We will therefore make the conservative 
approximation that T w = T D1 . 

Other 2^2 processes, for example those involving an Si in the initial or final state, 
have been neglected since their main contribution is at T > Mi where they would act to help 
create an initially thermal population of Si j2H] • These processes would tend to increase the 
leptogenesis efficiency in scenarios with zero initial abundance of Si, when K <C 1. The 
bounds derived later will depend on the leptogenesis efficiency at large values of K 3> 1, 
therefore we expect these processes to have a negligible impact. The same point can be 
made with regard to thermal effects, which would kinematically block the decays of Si 
at temperatures T > Mi [22]. As we consider scenarios in the large K regime, the baryon 
asymmetry is predominantly determined by processes at T < Mi, leading to only small finite 
temperature corrections to our T = estimates. 
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Figure 3: (Left panel) Leptogenesis efficiency, n defined in \29\i . versus K, for thermal and 
non-thermal initial abundance of S (S) and for various Br. (Right panel) Area in the M l7 
(bJh) n parameter space allowed by successful baryogenesis when (hj,h 1/ ) 11 = (h p hj,) and 
a = v = 175 GeV. 



Clearly, our treatment of the BEs is approximate. We have also neglected lepton flavour 
effects, coming from the charged lepton Yukawa couplings [2E] and the neutrino Yukawa 
couplings [21] • These effects could be large for scenarios with Mi < 10 12 GeV, when the 
charged lepton Yukawa couplings are in thermal equilibrium, or when there is only a mild 
hierarchy in the MjS. These effects could lead to differences in the final baryon asymmetry of 
up to an order of magnitude, however they are highly model dependent and would require a 
study beyond the scope of this paper. Since the purpose of our discussion here is to provide an 
existence proof and a reasonable estimate of the baryon asymmetry in very general scenarios 
our approach is expected to be accurate enough. 

After solving the BEs we can still parameterise the baryon asymmetry using the efficiency 
k defined in equation (|2*9j) . k will depend on the values of Mi, (h^h^n and (h p h_)n. Fig. [3] 
shows the dependence of k on K, for various initial conditions and various values of Br. In 
cases with no initial abundance of Si we see that the maximal efficiency k ~ 1 is indeed 
reached at K ~ 1. 

With either an initially thermal abundance of Si, or no initial Si, the behaviour of k for 
large K ^> 1 is the same, and for K > 20 is well fitted by the power-law 

Although much larger CP-asymmetries are required to produced the observed baryon asym- 
metry, this large K, or 'strong wash-out' regime clearly has the advantage of being insensitive 
to initial conditions. Fig. El also shows the behaviour of k for differing Br (or effectively the 
ratio h v : h p for small 5r). Br = 1.98 or Br = 0.02 corresponds to a factor of 10 difference 
between h u and h p . We see that as Br departs from 1 the efficiency for large K increases 
slightly. This effect is due to the less efficient wash-out via inverse-decays in these cases, as 
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can be seen in the BEs flSOjl where the second term in the equation for t)al responsible for 
the wash-out of the asymmetry via inverse decays is also dependent on Br. 

On the right hand side of Fig. 3 we show the area in the Mi, (h"'"h) 11 parameter space 
which is allowed by successful baryogenesis when (hj,h^) 11 = (h p h p ) n and a — v — 175 GeV. 
The points plotted correspond to numerical solutions of the BEs with the CP-asymmetry 
set to the maximum allowed by the bound (j2"Tj) . where the final lepton asymmetry would 
result in a baryon asymmetry equal to or exceeding the observed one. On the plot we also 
superimpose m iso-contours. 

If we take the most representative natural scenario, namely (h^h l/ ) 11 = (hph p ) n ~ 1 and 
the reasonable assumption that rh = 0.05 eV for hierarchical light neutrinos we can use the 
large K fit to the efficiency ()32j) and the bound on the CP-asymmetry (|2~T]) in conjunction 
with the value of the observed baryon asymmetry to set an approximate lower limit on er. 
We find that unless a > 0.1 GeV the L VR asymmetry produced is insufficient to explain the 
observed baryon asymmetry. Notice that this bound depends on several assumptions, in 
particular that the heavy Si are hierarchical in mass. Furthermore, the requirement that the 
universe reheats enough to thermally produce the S\ at the end of inflation leads to an upper 
bound on Mi of the order of Tr#, the reheating temperature. This leads to the approximate 
upper bound a < 2 TeV (T RH /10 ie GeV). 

In summary the scale of the spontaneous symmetry breaking of U(l)v is bounded, 

0.1 GeV < a < 2 TeV (^^y) , (33) 

and we can therefore conclude that an electroweak scale a is both natural and compatible 
with successful Dirac leptogenesis. 



4 The Higgs sector 
4-1 The potential 

The complete potential of the neutral scalar fields under consideration reads, 

V(H,<I>) = fi 2 H H^H + fil<t>^ + X H (H j H) 2 + X^^) 2 -r]H^H<t>^ , (34) 

where all the parameters are real. For the following, we denote H = H°. Notice that linear 
or trilinear terms do not appear thanks to the phantom U(1)d symmetry. After U(1)d 
is spontaneously broken, the field $ develops a vev a, which through the link 77-term in 
(134}) . forces the Higgs field H to also develop a vev v, triggering the "observed" electroweak 
SU(2)l x U(l)y symmetry breaking 2 . Expanding the fields around the minimum we obtain, 

H = v + -=(h + iG) , $ = a + -=(</> + U) . (35) 
v 2 v 2 



2 Notice that the limit — » is attainable and causes no problem for electroweak symmetry break- 
ing. However, we cannot justify a possible absence of the /j,^-term from (|34|l by using symmetry or other 
arguments. 
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While the Goldstone boson G is eaten by the gauge bosons, the same is not true for the 
remaining massless Goldstone boson J. Furthermore, the fields h and 0, under the influence 
of the 77-term, mix and become two physical massive Higgs fields Hi,i = 1,2 with masses 
m H2 > m H1 , 



and mixing angle 



tan 26 = - — . (37) 



The limits v ^> a and v <C a lead to the usual SM scenario with an isolated hidden sector. 
However, bear in mind that these limits require an un-naturally small 77 and will not be 
simultaneously compatible with neutrino masses and baryogenesis ()33)1 as advocated in the 
previous section. The case ^<1, with phenomenology resembling that of the SM, looks like 
an unnatural corner of the parameter space from the perspective of the Minimal Phantom 
SM. However, the exact limit 77 = is preserved under radiative corrections. The most 
interesting scenario is the most "natural" one when we require small Dirac neutrino masses 
and leptogenesis: tan# ~ 1. If the phantom sector is responsible for electroweak symmetry 
breaking then the natural choice of parameters, when taking into account the positivity 
constraint A# A$ > ?7 2 /4, is 

Xh ~ A$ ~ rj ~ 1 , tan# ~ 1 , tan/3 = v/a ~ 1 . (38) 

Under these conditions, succesful EW symmetry breaking happens only when /i| < 3 . This 
naturalness condition is supported by an enhanced symmetry $ <-» H of the potential (|3~IJ) 
which is broken however by the Yukawa couplings and the U{l)y- The Higgs potential of 
has been studied in [2H| in the context of Majoron models j3U] where a is an arbitrary 
vev. This study was mainly confined to LEP collider signatures. It is therefore interesting 
to update the phenomenology of this Higgs sector after the LEP era and in light of the 
forthcoming LHC experiments and the condition ()38j) . 

From (jnnj) we see that h = OnHi and therefore the couplings of Higgs bosons Hi to 
fermions and gauge bosons will be reduced by a factor On relative to their corresponding 



SM ones. It is almost obvious from ()34I35I36|) that Hi will couple to the "invisible" massless 
Goldstone pair J J. The situation is completely different to the SM where the H —* bb mode 
dominates for relatively light Higgs masses ^160 GeV. Here we find that in this mass range, 
the decay rate Hi — > J J relative to Hi — > bb reads as, 

^ = lf^i) 2 tan 2 /5tan^, (39) 

T(i7 1 ^bb) 48 V m b J K J 



^ - g = lf^) 2 tan^cot^ 
T(H 2 -> bb) 48 V m b J H 



(40) 



3 At this stage this has to be put in by hand i.e. we cannot offer a radiative mechanism to explain this. 
Dynamical EW symmetry breaking by fourth generation condensates has been discussed in |28| . 
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Figure 4: Dominant branching ratios of the two Higgs bosons decaying to invisible J J, and 
to SM particles bb, WW , ZZ and tt as a function of the Higgs boson mass H\ (left panel) 
and H2 (right panel) and 9 = 7r/4. In every case the mass of the partner Higgs boson is 
displayed in the upper horizontal axis for comparison. The shaded area is excluded by LEP. 

Therefore, from (J38I39I40|) we see that a light Higgs boson will decay dominantly to 
invisible J. J as long as it is heavier than 60 GeV. The existence of the massless Goldstone 
boson J diminishes the Higgs boson decay into a bb pair. LHC experimenters should be aware 
of this situation which arises in a very simple and natural extension of the SM! On the other 
hand, the Higgs decays to SM-vector bosons and fermions respectively, as 

Y{Hi^VV) = T{Hi^VV)\sMxOl , (41) 
T(H^ff) = r(Hi^ff)\ SM xO^, (42) 

with V = Z or W. Analogous formulae are valid for the cross sections a(e + e~ — > V* — > VHi) 
and o~(pp — > V* — > VHi). In the year 2001, the LEP collaboration presented bounds on an 
"invisible" Higgs boson mass Following this analysis, certain Higgs boson mass values 
are excluded function of the parameter 

? ^ a(e+e -> HZ) ^ g ^ invisible) = Ol x Br(# -> invisible) . (43) 
a(e+e HZ)\sm 

For ^ 2 = 1, LEP excludes Higgs boson masses up to its kinematical limit, m# < 114.4 GeV. 
This bound changes only slightly in our case. To make the above discussion more concrete 
we focus on a representative example with the natural choice of parameters given in ()38jl . 
i.e. 9 = (3 = 7r/4 and all couplings equal to one. In this case the production cross sections 
and branching ratios to gauge bosons, given by (|41j). amount to half of the SM predictions. 
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Thus £ 2 ~ 1/2 for Br(iJ — > J J) ~ 100% which is the case for Higgs masses in the region 
70 < rriHi ^ 160 GeV. This is nicely summarised in Fig. EJ where we plot the branching 
ratios for both Higgs bosons, Hi and H 2 , decaying into invisible Goldstone bosons J and to 
other SM-like particles. In general there are two more "invisible" decay "leaks" not depicted 
in Fig. EJ the first is H 2 — * HiH\ which is kinematically forbidden for our choice of 9 = 71 /A 
and the other is Hi — > vv which is proportional to the neutrino masses and is therefore 
negligible. It turns out that for mm ^ 160 GeV the LEP parameter £ 2 > 0.4 and therefore 
LEP [HI] excludes a light invisible Higgs boson with a mass mux ^ 110 GeV. This also sets a 
lower bound on the partner Higgs boson mass, mu2 ~ 191 GeV, which is now forced to decay 
only to visible particles WW, ZZ, tt. A search for the latter would follow the SM-type plan, 
looking for qqH 2 — > qqWW^ or gg — > H 2 — > WW^*\ZZ^*\ modes at the Tevatron and 
LHC. Regarding these channels, the only difference here is that production cross sections 
and decays are reduced by half. 

It is apparent from Fig. EJ that there is a mass region 

110<m m <160 GeV, (44) 

where Hi decays to invisible Goldstone bosons with Br(i/i — > J J) > 90%. The question is 
how can we identify this invisible Higgs boson at the LHC? This question has been studied 
extensively in the literature [321 EHE E3] • The purely invisible Higgs Hi can be searched for 
at the LHC through the Z + Hi and/or the W-boson fusion channels. Our analysis closely 
follows the results of [Bl| for the Z[— > l + l") +H mv production mode, where we multiply their 
S/\/B by a factor of 1/2 because of (j4*Tj) . We find that for an LHC integrated luminosity 
of 30 fb _1 the signal significance for the invisible Hi with a mass of 120 (140) [160] GeV 
is 4.9<t(3.6o") [2.7a] respectively. Although these results refer to the case where 8 = tt/4 the 
situation is rather generic in the region of ()38|). Although the above analysis for the Higgs 
boson decay to invisible is very important to identify the nature of the phantom sector one 
may also identify the light Higgs boson when mux ^ 140 GeV through the conventional 
Hi — > 77. We have not made a detailed study for this mode and we believe that it is 
worth further investigation. For mui > 160 GeV, Hi decays mainly to WW and ZZ since 
the Br(if — > J J) ~ 1/13 is suppressed. Notice however that mui ^ 200 GeV is rather 
disfavoured by the EW data as we will see shortly. 

In conclusion, the phantom sector of (p^fj) allows for Higgs decays into invisible Goldstone 
scalars and to visible gauge bosons. A situation may arise where LHC experimenters could 
detect Hi with mux 120 GeV through the invisible Z{-^ + Hi mode and H 2 with 

m H2 — 200 GeV through its production and decays in association with gauge bosons. 

4-2 The p-parameter and other observables: 

In a model with only Higgs doublets and singlets, the tree level value for the electroweak 
parameter, p = m^/Mf cos 2 9w, is automatically equal to one without further adjustment 
of the parameters of the theory. The correction to the parameter p, denoted by Ap, appears 
at one loop level. For the model at hand, the phantom singlet $ will affect gauge boson loops 
through its ^-mixing term with the observable Higgs field H. Then it is straightforward to 
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calculate the Higgs contribution to Ap It reads, 



Ap 



H 



3(jr p 

8y/2n'' 



rn 



In 



rn 



Hi 



m% In 



rn 



2 n 

Hi 



rn 



(45) 



where O is the orthogonal matrix in (}36|) . We can establish a useful connection between this 
formula and the SM one. Note that from the similarity condition of the rotation matrix, 
T m 2 = diag(m 2 Hl , m H2 ), with m being the 2x2 Higgs mass matrix, we read the following 
identity, 



m 2 m 0l = A\ H v' 



rn 



H ) 



(46) 



=1,2 



where m# is the SM Higgs boson mass expression. It is easy now to simplify our expression 
for Ap, by Taylor expanding ()45j) around m H , 



(47) 



where f(x) is a continuous function and f'(x) denotes its derivative with respect to m 2 Hi . 
Using (|46|) and the orthogonality condition T = 1, the second term in (|47p vanishes 
identically, leading to 



3G f 



8V% 



7T Z 



rn 



In 



in 



n 



rn 



9 , m M 

m\ In— f 



rn 



w 



(4? 



which is just the SM Higgs contribution to Ap. One arrives at the same conclusion for 
the S, T and, U parameters. Assuming that the Higgs contributions to the non-oblique 
corrections follow the same pattern, we can use the electroweak constraint on the SM Higgs 
boson mass, m# < 194 GeV at 95% CL jHEj in order to set constraints on the Higgs boson 
masses and mixing angle of this model. Thus, comparing (|45|) with (}4*5j) we arrive at 

(49) 



cos^logtm^) + sin 2 fllog(rr4 2 ) < log(194 2 GeV 2 ) (at 95% CL.) 



In the case of our working example 9 = 7r/4, this translates into mjjim^ 2 < 194 2 GeV 2 , e.g, 
m Hi ~ H5 GeV and ^5 327 GeV. These bounds have to be combined with the LEP 
bounds on the Higgs masses derived in the previous section. 

It is apparent that the inclusion of the phantom singlet field <3>, does not affect the GIM 
mechanism jHZj which is responsible for the absence of tree level flavour changing neutral 
currents (FCNC). One may think of Higgs mediated contributions to rare B-decays at loop- 
level. The most striking one would have been: B(or Y)-meson decays to invisible, B — > J J. 
Alas, the amplitude for this decay is proportional to £\ OnOi2 which vanishes because the 
matrix O is orthogonal. This is a kind of GIM mechanism suppression in the Higgs sector. 
Other Higgs mediated contributions to observables like B — > p + pT or to the muon anomalous 
magnetic moment, g — 2, will follow the SM prediction thanks to the relation (|47j) and the 
bound on the light Higgs mass (j44j) . In conclusion, the minimal singlet phantom sector does 
not change the FCNC predictions for processes existing in the SM. 
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5 Cosmological and astrophysical constraints 



Let us finally address the implications for cosmology and astrophysics. The presence of the 
massless Goldstone boson, J, has interesting consequences which need to be analysed. 

During the expansion of the Universe, a critical temperature is reached below which the 
U(1)d symmetry is spontaneously broken. As we have already explained in the previous 
section, the field $ then develops a non- vanishing vev and its real part, <j), mixes with the 
real part of the SM Higgs field, giving rise to two scalar Higgs mass eigenstates, Hi. On the 
other hand, J survives as the massless Goldstone boson. From (j34|35|36|) Hi couples to the 
Goldstone pair J J as 

(\/2C l 1 / 2 

-Cj D 1 F) tan /3 O l2 m 2 H . H t JJ . (50) 

The Goldstone bosons are then kept in equilibrium via reactions of the sort J J <-> //, 
mediated by Hi. However, since the amplitudes for these processes are suppressed due to 
a GIM-like mechanism which stems from the orthogonality condition, EjOaOi2 = 0, of the 
matrix O, J falls out of equilibrium before the QCD phase transition and remains as an 
extra relativistic species thereafter. 

The presence of relativistic particles, apart from neutrinos, is strongly constrained by 
Big Bang Nucleosynthesis (BBN), since they alter the predicted abundances for the light 
elements. Namely, additional relativistic particles would increase the expansion rate of the 
Universe, leading to a larger neutron-to-proton ratio and therefore a larger 4 He abundance. 
The allowed number of extra relativistic degrees of freedom is usually parameterised by 
the effective number of neutrino species, N e ff = 3 + AN U . Observations of the primordial 
4 He abundance, combined with the CMB determination of the baryon-to-photon ratio yield 
N u = 3.24 ±1.2 at 90% CL [23122], and a similar upper bound on N v can be derived from 
analysis of the CMB and large scale structure |H3] . Notice that this does not pose a problem 
for Js. Since they decoupled at a temperature above the QCD phase transition, when 
g* > 60, their temperature at BBN, Tj, is smaller than that of neutrinos and photons, T. 
Namely (Tj/T) 4 < (10.75/60) 4 / 3 . This is equivalent to an increase in the number of neutrino 
species of just AN U = 4/7 (Tj/T) 4 < 0.06, well in agreement with the above-mentioned 
constraints. 

On the other hand, J also induces the decay of heavy neutrinos into the lightest one, 
vh J vl- Their interaction is described by the following effective Lagrangian, 

C Jvv D —ivj 5 vJ, (51) 
a 

with m„ a diagonal matrix which contains the physical neutrino mass eigenstates. Since 
the decay does not include any photon in the final state, some potential cosmological prob- 
lems associated with radiative neutrino decays (e.g., contributions to the diffuse photon 
background and distortions to the cosmic microwave background black body spectrum) are 
avoided. However, limits on the non-radiative decay of neutrinos (usually expressed as upper 
bounds on the Jvv coupling) can be derived from solar neutrino observations (gj uj? % 10 -5 ) 
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[H], meson decays {g 2 j uP ^ 10 4 ) [12], as well as from preventing overcooling in supernovae 
(which exclude a range around g 2 JuU < 10~ 10 , although the bounds are model dependent) 



As we can read from (jSH), in the present model gj u p = m u /a. For natural (electroweak 
scale) values of a, and with m u < 1 eV we obtain g 2 JuD < 10~ 22 , thus fulfilling all the afore- 
mentioned constraints. 

Alternatively, a calculation of the heaviest 'light' neutrino lifetime, in the case of hierar- 
chical neutrino masses, yields 

— = ~ <?(10 13 ) s/eV, (52) 

where a ~ 100 GeV and m v ~ 0.05 eV, has been used. A similar value is obtained in the 
opposite limit, when neutrino masses are quasi-degenerate. The heaviest 'light' neutrino is 
therefore extremely long-lived and escapes all the constraints on neutrino decays. 

Finally, notice that the Goldstone boson J couples very weakly to electrons, through 
one loop diagrams, with strength g e j ~ Gi?m 2 m e /cr, and it does not affect the evolution 
of stars [13] if g e j < 10" 12 which implies that a > 5 x 10" 17 GeV. Also, the emission of a 
Goldstone pair mediated by (virtual) Higgses is negligible in our model, once more due to the 
orthogonality condition of (jHBj) . and the strong constraints on the Higgs couplings, obtained 
from studies on star evolution [13], are trivially fulfilled. 



6 Conclusions 

We have proposed the minimal, lepton number conserving, SU(3) c xSU(2)LxU(l)y gauge- 
singlet, or phantom extension to the Standard Model leading to naturally small Dirac masses 
for the neutrinos and baryogenesis through Dirac leptogenesis. The extension is natural in 
the sense that all couplings are either of 0(1) or strictly forbidden. 

Spontaneous breaking of a global, phantom sector U(l)£> symmetry will trigger elec- 
troweak symmetry breaking. The scale of this phantom sector symmetry breaking is con- 
strained to be around the electroweak scale by the simultaneous requirement of successful 
Dirac leptogenesis and small light neutrino masses. In this model, small Dirac neutrino 
masses arise through a mechanism very similar to the standard Majorana see-saw. The 
model can also be viewed as a very simple Froggatt-Nielsen scenario of the sort usually 
invoked to generate large hierarchies in the quark masses. 

Baryogenesis through Dirac leptogenesis occurs naturally in this model since the small 
effective Yukawa couplings of the left and right-handed neutrinos prevent the left and right 
neutrino asymmetries from equilibrating once they are created. The initial neutrino asym- 
metry is created via the out of thermal equilibrium decays of the heavy Dirac particles Si 
and Si, in analogy with conventional Majorana leptogenesis. 

A Davidson-Ibarra-like bound on the CP-asymmetry in the Si, Si decays exists when 
their masses are hierarchical. This bound, in conjunction with information on the efficiency 
of leptogenesis extracted from the solution of the Boltzmann equations allows us to place 



17 



a lower bound on the vev of the phantom sector, SM gauge-singlet $, such that the asym- 
metry created in Dirac leptogenesis is enough to explain the observed baryon asymmetry of 
the Universe. Assuming a hierarchical light neutrino spectrum and a hierarchical Si mass 
spectrum, we find that 

a > 0.1 GeV. (53) 

Making the further assumption that leptogenesis must proceed after the thermal production 
of the Si following a period of inflation leads us to an approximate upper bound on Mi and 
therefore a 

^ 2TeV (ife)' (54) 

Thus we find that an electroweak scale a is simultaneously compatible with both light neu- 
trino data and successful Dirac leptogenesis. Significantly, an electroweak scale a is also 
required by our naturalness criterion, since the mixing of the $ and the SM Higgs is ex- 
pected to be maximal. 

The addition of the phantom sector scalar which mixes with the SM Higgs, introduces 
an additional massive Higgs boson. After the breaking of the global U(l)o, we are also left 
with a massless Goldstone boson, J. This Goldstone boson couples to the two physical Higgs 
bosons introducing an additional, invisible decay mode for the Higgs Hi —> J J. This decay 
mode suppresses the branching ratio of Higgs to bb and instead both Higgs bosons decay 
dominantly to invisible J J and/or to vector boson pairs. We discuss in detail a natural 
scenario with a representative mixing angle (9 = vr/4) and estimate that with 30 fb" 1 of 
integrated luminosity the LHC could find the invisible Higgs with a mass 120 GeV, at a 
significance of 4.9 a. In addition, at the same time a significantly heavier, partner Higgs could 
be found with a mass 200 GeV through its vector boson decays. Interestingly, electroweak 
constraints suggest an upper limit of ~ 250 GeV for the mass of the Higgs bosons. 

The model passes relevant FCNC constraints thanks to a GIM-like mechanism. Cosmo- 
logical bounds on the number of relativistic species at BBN are also fulfilled, due to the 
Goldstone boson decoupling before the QCD phase transition. Finally, astrophysical con- 
straints on the Goldstone couplings from neutrino decays and stellar evolution are trivially 
satisfied. 
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